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O 

^ We study the prospects of controhing transmission of broadband and bi-chromatic 

^ laser pulses through turbid samples. The ability to focus transmitted broadband 

light is limited via both the scattering properties of the medium, and the technical 

^^ characteristics of the experimental set-up. There are two time scales, given by pulse 

c/^ stretching in the near- and far-field regions, which define the maximum bandwidth 



of a pulse amenable to focusing. In the geometric optics regime of wave propagation 

in the medium, a single set-up can be optimal for focusing light at frequencies uj 

and nuj simultaneously, providing the basis for the 1 + n coherent quantum control. 

^ Beyond the regime of geometric optics, we discuss a simple solution for the shaping. 
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which provides the figure of merit for one's ability to focus simultaneously several 



^ transmission modes. 

> 

O I. INTRODUCTION, 

in 

O This work is motivated by the goal of applying quantum coherent control techniques in 

^^ turbid samples. One of the basic ideas of quantum control involves focusing laser fields 

.^ of the frequencies uj and nuj onto an object with nonlinear response. Interference of the 

;-H excitation pathways due to each frequency component creates an asymmetric excitation 

[1]. This effect has been demonstrated in a number of theoretical and experimental works 

devoted to generating directed currents, control of absorbtion and propagation of light, 

and breakup processes in various physical and chemical systems [15]. Extensions of this 

principle, based on applying ultrafast laser pulses with controlled broad spectrum, have lead 

to numerous applications in control of quantum evolution, quantum information processing, 

spectral characterization, detection, microscopy and manipulations with microscopic and 



nano-scopic objects [1-4, 6-8]. We are interested in both the "1 + n" scenario and control 
with shaped ultrafast pulses. This task requires an ability to focus either bi-chromatic or 
broadband laser pulses with shaped spectrum in space in time. 

As a laser pulse is applied to a turbid sample - such as ground glass, biological tissue, 
paint, suspension, plastic, etc - its temporal and spatial structure breaks down [9-13]. In 
space, a coherent beam breaks into a multitude of speckles, so that spatial focusing is 
destroyed. In the spectral domain, the spectrum at each point in space can be strongly 
modified, so that the pulse shape is destroyed. The two effects are related, and each of them 
is deleterious for coherent control. 

This paper analyzes control of transmission of light with multiple frequency components 
in turbid samples, with the goal of designing quantum control experiments. For narrowband 
light, the corresponding technique [14] has recently lead to a breakthrough in focusing and 
manipulating laser beams in opaque samples [15-18]. The method is based on using a 
two-dimensional phase mask for the spatial correction of the wave front. We analyze the 
capabilities of this approach for spatio-temporal shaping of ultrafast laser pulses. While the 
first tests have demonstrated the great potential of the method for temporal focusing [19- 
21], efficiency of control over the broad bandwidth of ultrafast pulses needs to be thoroughly 
understood. Indeed, an experimental set-up optimized for controlling transmission at one 
given frequency, may not be suitable for another [22] . A set-up built to focus light at many 
frequencies simultaneously may be far from optimal for each individual spectral component. 
This work questions the fundamental limits of controlling broadband transmission through 
an opaque sample [15-19, 21, 22]. 

We find that the the ability to focus transmitted broadband light is limited via both the 
scattering properties of the medium, and the characteristics of the Spatial Light Modulator 
(SLM) used to modify the incident wave front. There are two time scales, given by pulse 
stretching in the near- and far-field regions (defined further in the text), which set the upper 
limit of the bandwidth of a pulse that can be focused. Their consideration suggests an 
optimization of the experimental set-up. In the geometric optics regime of wave propagation 
inside the sample, a single set-up can be optimal for focusing light at frequencies uj and nuj 
simultaneously, providing the basis for the 1 + n coherent quantum control, as demonstrated 
by our numerical simulations. Beyond geometric optics, i.e. when multiple interference can 
not be neglected, there is a simple figure of merit for one's ability to focus simultaneously 



several transmission modes in space. We also discuss a potential ability of using an opaque 
sample for shaping broadband spectrum, effectively replacing the dispersion element in the 
conventional pulse shaper [19, 21]. 

The rest of the paper is organized as follows. In the next Section we describe the implied 
experimental set-up, formulate our task in details, and describe the numerical simulations 
used throughout the text for illustration purposes. In Section III we neglect dispersion and 
backscattering, and solve the problem in the geometric optics regime, where the typical 
scatterer size is bigger than the laser wavelength. Thus we find the bounds on focusing 
imposed by the finite modulation depth of the SLM. In Section IV, we extend the description, 
including the effects of dispersion, finite spatial resolution of the phase masks's pixels, and 
of focusing of a laser pulse in time and at an angle. We also discuss focusing of broadband 
pulses in time vs. focusing in space. The general case, which goes beyond the geometric 
optics regime, is considered in Section V, where we discuss the scaling of the problem, and 
a simple strategy for using SLM to control simultaneously several independent transmission 
modes. In the last Section we summarize the findings of this paper. 

II. SET-UP. 

Fig. 1(a) shows the general set-up according to Refs. [15-19, 21]. The wavefront of a 
laser beam is modified by a two-dimensional SLM, whose pixels add a phase to the incident 
wavefront. The beam is then sent onto the scattering sample. Such a configuration allows 
the optimization of spatio-temporal focusing in either the near- or far field. 

Below, the "near field region", Eneari^^y^^)) corresponds to the output surface of this 
sample. The other, "far field" region, with the field distribution £", is at infinity along the 
z axis For a spatial harmonic transmitted at an angle ^, 

E{k^6) = / Eneari^^y) exp[i/cx siu 0] rfx rf^ (1) 

In this paper we concentrate on scattering that is sufficiently treated in the eikonal regime. 
We limit our consideration to focusing in the far-field, since it allows for easier modelling. 
This would be equivalent to optimizing transmission into a particular spatial harmonic of 
Enear) which cau bc then focused with a lens. As we explain below, the temporal structure 
of the pulse remains largely undisturbed in the considered regime and we primarily discuss 
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FIG. 1: (Color online) Implied experimental set-up. In panel (b), the SLM and imaging lens are 
replaced by the image of the SLM on the input surface of the sample. In addition, the direction of 
light propagation is reversed. 

focusing in space. 

Numerically, we solve the scalar wave equation for the electric field amplitude in the 
parabolic approximation [9]. The random medium is modelled by a set of planes. Each 
plane modifies the wave front as if the light was passing through a thin glass slide (refractive 
index n = 1.51) with randomly placed "impurities" characterized by a variation An in the 
refractive index. An example is shown in Fig. 2(a). Here the glass slide is taken to be 10- /xm 
thick, and round Gaussian-shaped impurities of the 1/e radius a = 30-/xm are characterized 
by An < 0.2. In the calculations, we place several such planes one after another, separating 
them by regions of empty space. 

Such modelling is inspired by experiments with diffusors based on random arrays of 
microlenses, ground glass, and all other opaque materials with relatively large (at least 
several microns in size) impurities [23-25]. A single slide in our modelling creates a far- 
field speckle pattern, but does not strongly modify the pulse spectrum. An array of slides, 
placed one after another, modifies both the spatial and temporal structure of a broadband 
pulse. Although our modelling misses the effects of de-polarization and backscattering, it 
allows one to understand some of the most important aspects of random propagation in the 
regime of low to moderate scattering angles (small backscattering). At the same time, the 
calculations are fast, allowing us to look at many frequency modes. Numerical propagation 
at each frequency consists of applying a coordinate-dependent phase to the wavefront at the 
location of each glass slide, followed by free propagation between the slides. The latter is 
made by making Fourier transform into the wave vector space and applying a fc-dependent 
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FIG. 2: (Color online). (a). Distribution of impurities in one of the four planes modelled in the 
calculation. The area shown is 770 x 770 fim. (b,c) Near-field intensity and phase distributions 
for a Aq = 800 nm beam after passing through four diffusing planes separated by regions of empty 
space. In panel (c), intensity is shown by brightness, and phase, between and 27r, is shown by 
color, (d) A 25-fs input pulse exits the system in the far field stretched to 100 fs. (e) Far-field 
speckle pattern for 800-nm light. 

phase to each spatial mode. For a femtosecond laser pulse sent into the sample, the temporal 
shape is obtained as a Fourier transform of the transmitted spectrum. 

Fig. 2(b) shows the near- field intensity of a 200-/xm-wide beam which has passed through 
a set of five planes with a = 30 /xm impurities. Adjacent planes are separated by 30 /xm of 
empty space. In the regime of geometric optics (a :» A), the speckle pattern is mainly due 
to multiple random lensing. Fig. 2(c) shows the same beam, stressing the phase at the exit 
from the last plane. The phase pattern is shown for the wavelength A = 800 nm. The speckle 
pattern at each frequency is almost the same, except for a frequency dependent phase which 
corresponds to a different time delay of the pulse arriving at different x, y points. The zeroth 
spatial harmonic is Eq = J Eneari^^v) dxdy. A 25-fs pulse sent to the system stretches in 






FIG. 3: (Color online) (a) Phase mask applied by ISLM to correct the wave front shown in Fig. 2(c). 
(b) Near-field intensity and phase distribution at 800 nm after the phase correction, (c) Corrected 
far-field intensity distribution. 

the far field to about 100 fs, as shown in Fig. 2(d). Fig. 2(e) shows the far-field speckle pattern 
for a A = 800 nm beam in the absence of the wave front compensation. 

Here we propose to image the SLM onto the input surface of a sample. Hereafter, we 
refer to the image of the SLM as ISLM. In this geometry, maximizing the transmission 
from the 0-th to the 0-th spatial mode (roo) is achieved simultaneously for forward- and 
for backward- propagating beams. Thus each pixel of the phase mask must add to the 
backward-propagating beam an x, y-dependent phase such as to make the wavefront phase 
as flat as possible (Fig. 1(b)). 

Within the arrangement of Fig. 1(b), we shall use the term "near field" for the field in 
the ISLM plane - even if it is placed at some distance from the actual border of the turbid 
sample. 

III. THE ROLE OF THE SLM'S PARAMETERS. 

ISLM can be thought of as a thin transparent plate with variable refractive index 
^S'lm(^,?/)- a wavefront passing through it acquires the phase 



(Pslm{x, y) = kaQ - ka{x, y) 



(2) 



where the optical pathways a{) — a{x^ y) are defined by nsLM{x^ y) = no + AnsLM{x^ y). Fig. 3 
shows focusing of the Aq = 800 nm wave shown in Fig. 2, with an ISLM of infinite spatial 
resolution, and a phase modulation depth of 27r (0 < a(x, y < Aq). Panel (a) shows the mask 



a{x^y) ranging from (black) to Aq (white). Panel (b) shows the corrected wavefront. Flat 
phase of the near-field wavefront ensures that the wave is almost perfectly focused in the 
far field, as seen in Panel (c). Variations in the near-field intensity somewhat decrease the 
focusing efficiency, adding a broad low-intensity pedestal, invisible at the scale of Fig. 3(c). 

We begin by neglecting dispersion and backscattering, and considering propagation in the 
eikonal regime. The latter corresponds to impurities in the sample being large, |Vn| <C fc, 
where k is the wave vector [10]. Experiments using commercially available diffusors, ground 
glass, random arrays of waveguides, etc, may fall under this case. In the eikonal regime, the 
wave is composed of trajectories - "rays". Each ray propagates in accord with the laws of 
geometric optics, and carries the phase kS^ where S is the optical path. The surface of equal 
phase at each point is orthogonal to the ray passing through this point; intensity variations 
are due to the varying density of the rays. 

Assume that the SLM's image has sufficient spatial resolution, and that the SLM is 
optimized to focus light with the wave vector ko = 27r/Ao. What happens with a wave 
characterized by k = ko + A/c? If the maximum ISLM's depth a was infinite, then the phase 
fiattening would work perfectly at each frequency. Indeed, by imaging the SLM mask on 
the surface of the sample we can effectively build a ffat slab out of the sample and ISLM: 
for each fc, 

k X [ao + So + AS{x^ y) — a(x, y)] = const (3) 

where S{x^ y) = Sq + AS{x^ y) is the optical path of a ray passing through the point (x, y) 
at the ISLM's plane. 

However, in reality the modulation depth a can cover only a few wavelengths. Assuming 
(^max = Ao, we have for the compensated wave front at /cq 

fco X [ao + 5^0 + AS{x, y) - a{x, y)] = 27rj^^ + ko{ao + Sq) (4) 

where jxy is an integer which can vary from one point (x, y) to another. This is the situation 
shown in Fig. 3. The term k^So is a constant phase which does not inffuence focusing. For 
a different wave vector, we have 

{ko+Ak)x[ao+So+AS{x,y)-a{x,y)] = {ko+Ak){ao+So)+27TJxy+Ak{AS{x,y)-a{x,y)) . 

(5) 
The phase compensation (5) will work for any k if the maximum modulation depth amax > 
AS for most pathways. If, on the other hand, AS' ^ oimax) the compensation will not work 



as soon as Ak{AS — a) c^ AkAS exceeds n for many points (x, y). Thus 
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According to the Huygens Fresnel principle, 

{AS)=CTfar, (7) 

where Tfar shows how much a short pulse sent to the system is stretched in the 0-th spatial 
mode or, equivalent ly, in the far-field focus. Another way to see this fact is as follows. 
Consider two points, A and B, at the exit from the sample, such that Sa = Sb + {^S). 
When CW light of the frequency uj is sent into the system, the phase of the field at the 
points A and B differs by (/>ab(^) = {AS)(jj/c. At a different frequency, uj + 6(jJ^ 

(I)ab{oj + Suj) = (I)ab{oj) + Suj{AS)/c (8) 

According to Eq.(l), the complex values of the field from all near-field points are summed 
to produce a far-field speckle. One can see from Eq.(8) that the detuning 6(jJ = tic/ {AS) 
corresponds to the speckle pattern being significantly different from that at the frequency 
(jj. At this value of the detuning, constructive interference between the fields coming into 
the far field region from the points A and B turns into a destructive one, and vice versa. 
Therefore, the frequency correlation length of the far field speckle pattern is, approximately 
2Aujfar = 27Tc/{AS). This means that the transmitted spectrum in the far field consists of 
independent bands of the width 2A(jjfar- Equivalently, a very short laser pulse sent into the 
system stretches in the far field to Tfar = Tr/Aujfar- 

Comparing Equations (6) and (7), we see that compensation can only work for detunings 
A(jU = (jU — CjUq such that 

\AUJ\ <^^= Aujfar . (9) 

^far 

Figure 4 assumes the compensation mask shown in Fig. 3(a) applied to the sample discussed 
in Figs. 2 and 3. Panel (a) shows the intensity Icenter in the far-field focus in dependance on 
the field wave vector, calculated for a single realization of the random sample. As k is detuned 
from ko = 12500 cm~^, the focusing vanishes. The width 2Ak c^ 500 cm~^corresponds, up 
to a numerical factor of c^ 1.5, to Tfar = 100 fs, i.e. the pulse stretching seen in Fig. 2(d). 

Unexpectedly, Icenter in Fig. 4 increases again in the vicinity of k = 2ko. The effect is 
explained in the following way. If the condition (4) is fulfilled, then 

2ko X [ao + So + AS{x, y) - a{x, y)] = Anj^y + 2A;o(S'o + S^^slm) , (10) 
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FIG. 4: (Color online) (a) Laser intensity in the zeroth far-field mode in dependance on frequency, 
with the ISLM mask tuned for cjq. (b,c) Near-field wave fronts dX uj — I.Scjq and uj — 2uj{). 

and the phase compensation at twice the main wave number is again complete, as shown in 
Fig. 4(c). We see that in the simplified model — negligible dispersion, ISLM can spatially 
resolve the phase front, — an ability to focus bi-chromatic fields, and to perform "1+n" 
quantum control comes at no expense. An experimental set-up optimized to focus a laser 
field at frequency uj will also focus field at frequency nuj. 

Note that the peak amplitude at 2/co in Fig. 4 is slightly smaller than that at /cq. Indeed, 
the assumption that the phase mask is able to resolve individual pathways becomes invalid 
at the near- field caustics, where several rays intersect at the same point. This situation is 
mathematically similar to that of an SLM with limited spatial resolution, discussed in the 
next Section. 

Moreover, similar to the case of /c = 2/co, at A; = L5/co the phase of the compensated 
wave can only have two values, and tt, as seen in Fig. 4(b). Each part of the near-field 
wave front ~ that with the zero phase, and the phase equal to tt - yields a strong focus in 
far-field. The two foci interfere destructively. However, because of the random amplitudes, 
the destructive interference is not complete, and focusing at A; = L5/co is still better than 
that at the adjacent values of k. Reminiscent of fractional quantum wave packet revivals 
[26], such incomplete focusing happens at any k = P/Q x /cq with integer P, Q. 
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IV. ADDITIONAL BOUNDS. 

The above consideration remains valid if the goal is to optimize transmission into a spatial 
harmonic propagating at an angle 6 (or, equivalently, off-axis far field focusing). In this case, 
Eq.(4) turns into 

/co X [ao + 5^0 + AS'(x,?/) - a{x,y)\ = 27r j^y + ko{ao + So) + k^x sin (9 (11) 

where x is the coordinate in the ISLM plane. In a compete analogy with Eq.(lO), light 
with the wave vector 2/co will also be focused. In our numerical simulations, the spectral 
bandwidth Ak of the spatially focused light did not depend on 9. 

The ability of the scheme to focus several frequencies simultaneously depends on the 
sample's dispersion. Indeed, the above consideration is based on the assumption that light 
at each frequency propagates along the same set of rays. In another series of calculations 
we included the effect of dispersion, assuming that the samples are made of BK-7 glass [27]. 
We found that the focusing survives in the presence of dispersion: In our calculations, the 
focused intensity at the frequency 2(jUo decreases only by a factor of, approximately, 2-3. This 
number is small compared to the ^ 10^-fold increase in the intensity at the focus observed 
in the case of complete phase compensation. 

Finite size of the ISLM's pixels in the (x^y) plane does bring an important additional 

bound on one's ability to focus broadband light. If the ISLM grid can not resolve the phase 

variations in the incident wave front, then each pixel will be used to compensate the phase 

of the average field 

Eaviko) = J2 i^(A5,)e^^°(^°+^^-) (12) 

s 

where P{ASp) is the probability distribution for the pathways characterized by the length 
So + ASp averaged by a single pixel. Coarse graining over ISLM's pixel size limits the 
compensation fidelity. Suppose that a pixel is tuned to compensate the phase of Eav at the 
given position at the frequency koc. At a different frequency we have 

Eav{k) = Y^ p(A5'^)e^[^^o+^oA^,+AfcA^,] ^^3^ 

s 
The values of Eav correspondent to k and k + Ak differ drastically if AkASp ^ n for many 
pathways passing through the particular pixel. Thus the phase compensation will not work 
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FIG. 5: (Color online) Focusing as a function of frequency in the limit of low spatial dispersion, 
for detunings Auj exceeding 
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(14) 



where Tnear = c {ASp) describes stretching of the pulse in the near field, averaged over an 
area of the ISLM's pixel. 

The bounds on the focusing imposed by the SLM are summarized in Fig. 5. There can 
be several peaks of the focused field, each of the width 2A(jjfar^ under the envelope of the 

width AuJnear- 

Note that Aujnear is due to the pulse stretching at the ISLM's pixel. If one moves the 
ISLM plane away from the surface of the sample, then Auj^ear approaches Aujfar- Same 
happens if the pixel size is increased, or if interference of multiple pathways at each point of 
the ISLM plane becomes too strong. We numerically verified that the peak at 2cl;o in Fig. 5 
disappears if the ISLM pixels become so large that they can not resolve the phase variations 
in the scattered wave. 

The peak also disappears if the interference in the scattering process can not be neglected. 
In our calculations, this was achieved by reducing the typical size of the impurities while 
simultaneously increasing their number. This led to both higher scattering angles, and 
deviations from the eikonal regime which allows interpreting propagation of light via an 
ensemble of rays. Surprisingly, however, predictions based on the eikonal optics approach 
hold even for rather strong scattering: For light that has passed through 30 3-/xm thick 
planes with a = 2 /xm impurities, the phase compensation for the Aq = 800 nm far-field 
focusing still provided a noticeable focal spot at A = Ao/2 = 400 nm. 

Our consideration above refers to focusing laser pulses in space, but not in time. Note, 
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however that the abihty to spatiaUy focus broadband hght is bound by Aujfar - the spectral 
bandwidth of a pulse which is not strongly distorted in the far field. Thus we show that the 
above approach is limited to a spectral band where the temporal structure is not destroyed, 
and temporal focusing is not required. If the far-field stretching is substantial, one needs 
to assign different pixels to different bands, as described in the next Section. Only at that 
stage the question of temporal shaping and focusing adjusting the relative phases of several 
independent frequency bands - arises. 

V. THE GENERAL CASE 

In the general situation, many interfering pathways may lead to the same near-field point. 
As before, let Tnear characterize the stretching of the pulse at the exit from the sample, and 
^far ~ at infinity. Similar to the previous Section, the near-field speckle pattern changes 
at detunings exceeding Acj^ = n/ Tnear ^ and finite depth of the SLM's phase modulation 
prevents focusing at detunings exceeding Auf = n/rfar except for the frequencies cj^ related 
to (jUq by 

/ kjndS{kjn) = m X kodS{ko) . (15) 

As shown in Section IV, if Tnear < ^far^ it makes sense to place the ISLM in the near field 
with respect to the random sample. Note that such situations include those where the pulse 
is significantly modified after passing through the random sample, both in space and in time. 

In the case of stronger scattering, when the geometric optics -based model is inapplicable, 
one must view the far-field focusing as phasing together random phasors corresponding to 
different scattering channels [14]. This regime is mathematically similar to that of large 
ISLM pixels, discussed in the previous Section (Eqs. (12-14)). Below we briefiy discuss what 
scaling should be expected for focusing broadband pulses in this case. 

Let us assume that the far-field transmission spectrum within the bandwidth of the laser 
pulse consists of M un-correlated bands of the width Acj/, and that the laser beam covers A^ 
pixels of the phase mask in the scheme of Fig. 1(a). In order to obtain figures of merit for the 
focusing capability, we assign N/M pixels of the SLM to each frequency band, in the way 
that is discussed below. Assuming a circular gaussian distribution for the field amplitude 
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after the sample [11, 14], 



P{ERe, Elm) = 7^—r exp 



^Re + ^Im 



(16) 



where P is the probability density, and Ejie, Ejrn are the real and imaginary field amplitudes, 
the intensity of the focused field at frequency u is [11] 



TT /N^^ 



In order for the focused spectrum to be controllable, this value must exceed the background 
due to the other pixels assigned to all other frequencies. The latter is obtained with the 
help of Eq.(16) as 

{hack{c^)) = {M-l)-l0. (18) 

Enhancement in spectral intensity due to the focusing is then [14] 

Once control over each spectral band of the width Aujfar is achieved, one can tune the 
overall phase of the field in each band by applying an extra phase to each phase mask's 
pixels controlling the mode. Through these phases, the spatially focused pulse can either 
be focused in time or be given any temporal shape allowed by the frequency resolution of 
A(jj and the number of pixels in the phase mask. Thus the system makes an analog of a 
conventional pulse shaper [28], with the dispersion element being replaced by the random 
sample [19, 21]. 

If the M spectral components are given equal phase, together they form a pulse that is 
M times shorter in time than each of the M components. Its intensity is M times higher 
than that of the incoherent sum of the components. Thus the maximum achievable intensity 
is 

times stronger than that of un-compensated light. 

VI. SUMMARY 

Coherent control of physical and chemical processes in turbid media requires availability 
of focused laser pulses with tunable temporal/spectral shapes. This, in turn, sets the task 
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of coherently controlling propagation of bi-chromatic and broadband laser pulses through 
turbid media. 

Recent works have shown that this task can be carried out by using phase masks to 
adjust the phases of different transmission modes. Each mode, centered at its own central 
frequency and having its own speckle pattern, can coherently contribute to the output 
field. By controlling the interference between the modes one can achieve the desired spatio- 
temporal focusing. In this sense, the experimental scheme shown in Fig. 1 is an analog 
of a conventional pulse shaper, with the dispersive element replaced by the turbid sample. 
Resolution of this turbid pulse shaper is set by the transmission properties of the sample 
[19], together with one's technical ability to control relative phases of the modes. The 
ability to shape pulses simultaneously in space and time, and to work with very narrow- 
band transmission modes, can bring new dimensions into experiments on coherent control. 

Most present-day experiments do not assume correlations between the phase patterns of 
different frequency bands, and work in the regime where the phase mask can not resolve 
phase variations within a single speckle pattern. In this case one can obtain the figure of 
merit for the efficiency of the spatio-temporal focusing of light by assigning a fraction of the 
phase mask to each of the independent frequency bands. This is done in Section V of our 
paper (Eqs. (19,20)). For M independent frequency bands, this leads to focused intensity at 
a single frequency scaling as tj^ oc 1/M^. If the phases of the frequency bands are set such 
as to produce a short pulse in the focus, its intensity scales as rit oc 1/M. 

An interesting regime arises in the case of moderately strong scattering and relatively 
large-size (above 2 jjlvh in our simulations with 800-nm light) scatterers. In this case, in 
agreement with the geometric optics approach, optimization of spatial focusing at frequency 
cjo automatically ensures that focusing at the frequency uj = nuj^ is also achieved. In this 
situation, " 1+n" coherent control must be available at no extra cost provided the relative 
phase between the two fields can be maintained. In addition, interesting phase structures 
arising at frequencies that are rational fractions of cjq call for further investigation. 

In the geometric optics regime, the efficiency of the spatial focusing is bound by the two 
time scales. A single set-up of the phase mask can only optimize spatial focusing within 
a single frequency transmission band, with the width given by 2A(jjfar = "^Tr/rfar (Eq.9), 
where Tfar corresponds to the stretching of an ultrashort pulse in the far field at the output. 
At the same time, there is an overall envelope of the focusing efficiency (Fig. 5). Its width 
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is given by Acj^ear = ^ linear- Here Tnear IS the duration of a pulse covered by the area of a 
single pixel of the phase mask in the geometry of Fig. 1(b), and /S^uOnear is the band with of 
the transmission matrix taken at one pixel. The separation of the two time scales suggests 
that one should choose the experimental set-up with the shortest Tnear- To minimize the 
pulse stretching in the plane of the phase mask, we proposed to image the SLM onto the 
input surface of the turbid sample. 

The intuition inspired by geometric optics remains valid if one considers far-field focusing 
at an angle, or if moderate dispersion of the sample is taken into account. However, using 
the same phase mask to focus at frequencies uj and nuo simultaneously becomes difficult 
if the pixels of the phase mask can not resolve individual near-field speckles. In this case 
/\ujnear approachcs AijUfar^ and pulse stretching in the near- and far field is the same. Then 
the geometry can not be optimized by placing the phase mask at any particular distance 
from the sample. Control over focusing of multiple frequencies can be achieved by assigning 
subsets of the mask to different frequency bands. 
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